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Abstract — Delay-Tolerant Networks (DTNs) have emerged as 
an exciting research area with a number of useful applica- 
tions. Most of these applications would benefit greatly by a 
reduction in the message delivery delay experienced in the 
network. The delay performance of DTNs is adversely affected 
by contention, especially severe in the presence of higher traffic 
rates and node densities. Many-to-Many (M2M) communication 
can handle this contention much better than traditional one- 
to-one communication employing CSMA. In this paper, for the 
first time, we analytically model the expected delivery delay of a 
DTN employing epidemic routing and M2M communication. The 
accuracy of our model is demonstrated by matching the analytical 
results against those from simulations. We also show using 
simulations that M2M communication significantly improves the 
delay performance (with respect to one-to-one CSMA) for high- 
contention scenarios. We believe our work will enable the effective 
application of M2M communication to reduce delivery delays in 
DTNs. 

I. Introduction 

Delay-Tolerant Networks (DTNs) JT], J5] are challenged 
networks where end-to-end connectivity does not always exist 
between nodes. One reason this may occur is that links get 
disrupted either due to interference or due to nodes moving 
out of range of each other. DTNs can be put to a variety of 
applications such as bringing Internet connectivity to villages, 
and reducing the pressure on cellular bandwidth by offloading 
delivery through WiFi and mobiles of people moving in the 
vicinity. Since an end-to-end path seldom exists between 
source and destination, most existing Internet protocols fail 
and instead, the store-carry-forward [1| protocol is used to 
transfer messages from source to destination through a num- 
ber of intermediate relays. The lack of an end-to-end path 
means message delivery delays are large. Hence, lowering the 
message delivery delays would be critical in enhancing the 
performance of most of these applications. 

It has been shown that contention adversely affects the 
delay performance of DTNs Q. This is because contention 
implies that transfer of a message is not always possible 
when one node meets another, owing to two of the following 
reasons: (i) the message contends with other messages in the 
buffer of the node to be chosen for transmission, and (ii) the 
node contends with other nodes in the vicinity for opportuni- 
ties to communicate. 



While contention has been recognized as an important 
factor that increases delivery delays, the problem of fighting 
contention and thereby reducing delivery delays in DTNs has 
not been adequately considered. For a lot of DTN application 
scenarios, it is possible for a number of nodes to come in 
the vicinity of each other. For example, people in a crowded 
public square, or villagers near a road where a vehicle is 
passing by. In these scenarios, as a result of contention, most 
nodes would not be able to transmit the messages they carry 
since existing access schemes such as Carrier Sense Multiple 
Access (CSMA) are one-to-one communication schemes and 
would allow only one pair of nodes within the region of 
contention to communicate at any instant of time. 

Many-to-Many (M2M) communication QJ, has been 
proposed to improve the capacity of Mobile Ad hoc NET- 
works (MANETs) and enables each node in a bounded group 
of nodes to simultaneously communicate with each other node 
in the group. Hence when a number of nodes come in the 
vicinity of each other, this precious contact duration is used 
much better, with more nodes getting a chance to communi- 
cate. Hence, M2M communication can be used to battle the 
effects of contention and therefore reduce message delivery de- 
lays. M2M communication has become feasible since the state- 
of-the-art technology enables support for both access schemes 
(such as Code Division Multiple Access (CDMA)) and loca- 
tion services such as Global Positioning System (GPS)) on a 
single IC chip 0]. 

In this paper, we analytically model the expected deliv- 
ery delay of a DTN employing epidemic routing and M2M 
communication. Despite the fact that the analysis has been 
done for epidemic routing, the framework we develop is easily 
extended to other routing protocols. The accuracy of our 
model is demonstrated by matching the results from analysis 
against those from simulations. Also, simulations show that for 
high-contention scenarios, M2M communication significantly 
outperforms one-to-one communication employing CSMA in 
terms of delivery delay. 

II. Related Work 

Initial analytical models developed for DTN routing per- 
formance study (6), 0, J8) worked under the assumption 
that whenever two nodes are in contact with each other, all 



messages could always be successfully transferred from one 
node to the other (i.e., they assumed both buffer capacity 
and bandwidth to be infinite). While papers such as [9] have 
modeled DTN performance with bounded buffer capacity, they 
have assumed that infinite bandwidth is available and hence 
that there is no contention. The motivation for not considering 
contention has been that DTNs are sparse networks and such 
sparsity yields negligible contention. However, this conjecture 
has been disproved with the help of simulations in works 
such as Q, lHJ. The authors in (3) show via simulations 
that, irrespective of whether the network is sparse or dense, 
the contention is substantial for high traffic rates; and also 
that the contention increases with an increase in network 
density. Realizing the importance of contention in the routing 
performance, the authors have attempted to include contention 
in the analysis of routing [3|, |10|. Their analysis assumes the 
one-to-one CSMA communication scheme. 

In our paper, however, we propose that in scenarios where 
there are areas of high node density (i.e., a number of nodes 
in the vicinity of each other), nodes must co-operate using an 
M2M communication scheme such as the one proposed in [4| 
and not compete as assumed by [3|. This would mean that 
instead of just one pair of nodes in such an area of higher 
node density communicating, we can have a larger set of 
nodes communicating with each other at the same time. This 
would significantly increase the probability that two nodes 
which come into contact with each other get to exchange 
messages, and thereby reduce the routing delay. The beneficial 
effect of using such an M2M communication scheme would 
be amplified in the presence of high traffic loads. 

III. Analytical Model 

In this section, we develop an analytical model for the 
expected delivery delay of a DTN employing M2M commu- 
nication. We follow the delivery of a particular message from 
source node to the destination node, and model the evolution of 
the state of the system with respect to this message. We first 
describe the Medium Access Control (MAC) protocol used 
for achieving M2M communication which logically divides 
the network into square cells (or) tiles. We then derive the 
mobility statistics for the Random Direction (RD) mobility 
model, where the network is logically divided into square tiles 
as described in the MAC protocol. Using these statistics, we 
model the evolution of the network as a stochastic process and 
thereby derive the expected delivery delay. 

A. MAC Protocol for M2M Communication 

We adopt the scheme proposed by Moraes et. al., in 10] for 
M2M communication in MANETs. While they have proposed 
this scheme to increase the capacity of MANETs, we adopt 
the scheme to reduce delivery delays in DTNs. For the sake of 
completeness, we briefly describe the scheme in this section. 
The network is logically divided into a number of square 
tiles. Each node, with the help of GPS, is aware of its 
current location (thereby the current tile it is in). GPS also 
helps in time synchronization among the nodes. Therefore, 



the communication in each tile is synchronized. In each 
tile, communication is performed in the form of sessions 
that have two phases. The first phase, called the discovery 
phase, is divided into slots. Each node in the grid that wants 
to communicate randomly picks one slot and broadcasts its 
presence to the other nodes in that slot, along with some 
other control plane information. The first bunch of nodes (say 
a-nodes) that successfully do so enter the second phase, 
which is the communication phase. The parameter a is chosen 
according to different factors such as the allowable receiver 
complexity. We assume there are sufficient slots to make the 
probability of collision negligible thereby ensuring that if there 
are lesser than a nodes in a particular cell, then all of them 
are chosen for communication with a high probability. Also, 
this means that when there are at least a nodes in a particular 
cell, the MAC scheme effectively chooses a random subset of 
a-nodes from the set of nodes in that cell. Each node that 
enters the communication phase is assigned a code on which 
it transmits data and a frequency on which it receives data. 
Using narrowband FDMA/CDMA, all a nodes are able to 
communicate with each other simultaneously (i.e., each node 
participating in the communication session can send up to 
bsw messages to every other node participating in the session, 
where bsw depends on the bandwidth available for M2M 
communication). 

B. Mobility Statistics 

1) Intermeeting Time (IMT): The IMT of any two nodes is 
defined as the time elapsed between two successive meetings 
of the nodes. Since we logically tile the terrain to make the 
MAC easier, the IMT is expected to increase since some 
communication opportunities are lost as a result of tiling 
the terrain. Hence existing derivations ifTD do not hold as 
such. We adopt the methodology used in |12| to derive the 
distribution of the IMT, show it to be exponential and derive 
its mean. We define the area covered by a node while moving 
along a path as the total area that came under the node's 
coverage when it moved along that path. It must be noted 
that a node's coverage area at any point is the cell it is in at 
that point. 
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Figure 1: Node performing RD motion 
Lemma 1: The average area A covered during an epochQ 

1 Under RD mobility, a node chooses its mobility parameters such as speed 
and angle of direction, and travels with the chosen mobility parameters for a 
duration of time drawn from an exponential distribution (this duration of time 
is called an epoch). After that, the node halts for a random amount of time, 
following which it chooses a new set of values for its mobility parameters 
and the process repeats. 
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under RD mobility over a square-tiled terrain is given by 
A = 2a 2 + where a is the length of a side of a cell, 
and L is the average distance covered in an epoch of the RD 
model. 

Proof: Consider a node starting an epoch within a cell 
at the coordinate (x, y) as illustrated in Fig Q] Since the RD 
model has a uniform stationary distribution, these coordinates 
x and y are drawn from a uniform distribution over [0, a]. Also, 
according to the RD model, the node chooses an angle 9 drawn 
from a uniform distribution over [0, 2tt], velocity v drawn from 
a uniform distribution over [V m i n , V max ], duration of epoch t 
drawn from an exponential distribution with average 4, and 
a halting time at the end of the epoch drawn from a uniform 
distribution with average T s . Let the number of vertical walls 
and horizontal walls crossed during this epoch be T) v and rjh, 
respectively. We have for 9 £ [0, ?], 



Hence we have 



and 
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We now derive the area covered in the epoch, by counting 
the number of cells that are traversed in that epoch. For most 
starting configurations, the node does not pass through corner- 
points of cells during the epoch. Hence, it is reasonable to 
assume that every wall crossed during the epoch is intersected 
somewhere along the side and therefore that each intersection 
represents crossover of the node into a new cell. Hence the 
number of cells n covered in the epoch for 9 £ [0, ^1 is given 

by, 



n = l + r/ v + rj h 
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Because of symmetry, the average number of cells covered 
in an epoch for 9 £ [0, ?] is the same as the average for 
9 £ [0, 27r]. Hence, the average number of cells n covered in 
an epoch is given by Eq. [T] which on solving gives, 



A = no 2 = 2a 2 
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Having estimated the average area covered during each 
epoch, we can proceed as in lfT2l to show that the expected 
hitting time0 Th is as follows: 



T h = 



N 



2a 2 + ^ 



L 



(3) 



In a similar way, as done in ifTZl . we can show that the meeting 
time T m is exponentially distributed with an average given by 
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(4) 



where v r d is a constant (1.27) for the RD model such that the 
expected relative speed between two nodes u^! is given by 



Vrel = VrdV, 



(5) 



and p m is the probability that a node is moving at any point 
in time given by, 

Finally we can show that the intermeeting time is also ap- 
proximately exponentially distributed such that its average Ti m 
equals T m . 

2) Contact Time ( CT): The contact time between two nodes 
is defined as the time elapsed between when the two nodes 
come within communication range of each other and when 
they go out of the communication range of each other. In this 
section, we derive the expression for the expected contact time 
(T con ). We assume that the probability of contact between two 
nodes lasting over multiple cells is negligible. When a contact 
occurs, two cases arise: either one node could be stationary or 
both nodes could be moving. Let the expected contact time in 
the first case be Tjf on and in the second case be 

Lemma 2: 



(1 -p m )T c s on + 



Proof: The probability that both nodes are moving given 
that one node is moving is p m and the probability that one 
node is stationary given that one node is moving is 1 — p m . 
The result follows from this. ■ 



(2) 



2 The hitting time is the time taken for a node starting from the stationary 
distribution to hit a randomly chosen static target. 



Lemma 3: 



where 



Ka 2 



log 



Vrr, 



log (\/2 + l) +(V2-1) 



and 



J ow = 2Ka 2 log (^-) \log (y/2 + l) + (1 - y/2) 



for if given by 



if 



Proof: When one node is stationary, the contact begins 
when the other node enters the cell and ends when the other 
node leaves (our assumption that contacts last only over a cell 
is especially valid in this case). Hence the expected contact 
time in this case T* on is given by the expected time elapsed 
between the moving node entering and leaving the cell. We 
assume the moving node enters through any particular wall 
of the cell (since the problem is symmetric w.r.t the different 
walls of the cell) with a uniform distribution on point of entry, 
velocity, and angle of entry. Since epoch lengths are assumed 
large, it is assumed that with high probability, the movement 
of the node through the cell is part of one epoch and hence a 
straight line. The node can leave either through one of the two 
adjacent walls or the opposite wall. I aw is the product of the 
expected time taken to leave through a particular adjacent wall 
with the probability that the node leaves through that adjacent 
wall. Similarly, I ow is the product of the expected time taken 
to leave through the opposite wall with the probability that the 
node leaves through the opposite wall. Hence, 



Tr n n law "i" law ~t~ low 2<I aw + I Q 



(7) 



I aw and I ow can be derived using geometric arguments 
which we do not present here due to space constraints. ■ 

Now we approximate the contact time distribution with one 
node stationary as an exponential distribution with expected 
value T£ on . Hence the contact time for both nodes moving is 
given by the minimum of the times spent by each node in 
the ce ll and is hence approximately exponential with expected 
value -f^- Therefore 



IV. Modeling Contention under M2M 
Communication 

In this section we develop a framework for modeling 
contention in a DTN employing M2M communication. We 
first derive the distribution for the number of neighbors of 
any node and then derive the expected number of messages in 
the buffer of a node. These results are then used to develop 
the contention model. 



A. Distribution of Number of Neighbors 

We have established that the IMT and CT are approximately 
exponential with means T im and T con , respectively. We hence 
model the process of nodes arriving within communication 
range of a particular node (i.e., entering the same cell as 
the node), and moving out of communication range of the 
node as a discrete time Markov chain, where the states are 
given by the number of neighbors at the beginning of each 
communication session. We assume that the probability of 
having more than N max neighbors is negligible. We also 
assume that since the duration of a communication session (r) 
is small, the probability of more than one arrival or departure 
in the duration of a communication session is negligible. 
Hence the transition probability P(n + l\n) of the Markov 
chain is the probability of one arrival and no departure and is 
given by, 



P(n + l\n) 
where 

and 
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In a similar way, P(n — l\n) is given by 

\eff T 



P(n+l\n) 



efS, 



eft. 



1! Ml! 
and P(n\n) is given by, 

P(n\n) = 1 - P(n + l\n) - P(n - l|n). 

The steady state distribution of this Markov chain gives us 
the probability distribution of neighbors. We assume in the 
subsequent sections that the probability that the number of 
neighbors N n hb equals j is denoted by vj. 

B. Expected Number of Messages in Buffer 

In order to derive the expression for the expected number 
of messages in buffer (E[B]), we note that by Little's Law, 
we have, 

E[B] = XE[W] (8) 

where A is the total traffic arrival rate at each node (assumed 
Poisson) and E[W] is the expected waiting time for a message 
in a node's buffer. 

We first derive the expression for the arrival rate A. We then 
derive E[W] in the next section. A is given by 

A Xgen A re / 

where X gen is the rate of generation of messages at that node 
and X re i is the rate at which relay messages (messages for 
which this node is going to act as relay) arrive at the node. 
To calculate A re ;, we assume that any given message spreads 
to nearly the entire network before being removed by the 



recovery mechanism. This is valid for the direct recovery 
mechanism we assume is used, which makes minimal effort at 
recovery. Hence, we assume that any message in the network 
is replicated at a node with a high probability before all its 
copies are erased from the network. Now, working under this 
assumption, we treat the network as a queue with messages 
generated at other nodes considered as arrivals (with an aver- 
age rate X gen {n to t — 1)), and messages delivered to a particular 
node under consideration as service completions (with an 
average rate \ re i). Clearly, for the number of messages in 
the network to be stable, the service rate A re ; must be at least 
Xgen{ntot — 1))- We approximate the value of A re z by this 
lower bound. 

1 ) Expected Waiting Time for Epidemic Routing: The ex- 
pected waiting time depends on several factors such as the 
chosen mobility model, routing protocol, and the recovery 
mechanism. In order to derive the expected waiting time for 
epidemic routing, we assume direct recovery is used, i.e., a 
node will only drop a message after meeting the destination 
and either transmitting the message to the destination itself 
or learning from the destination node that the destination 
node has already received the message. We also assume that 
the scheduler at each node works such that whenever the 
destination of some messages is encountered, those messages 
are given top priority. We also assume that if contact is 
established with the destination of some messages in the 
buffer, the bandwidth available is sufficient to transmit all the 
messages to the destination. This is valid since we assume 
symmetric traffic generation, i.e, for any source node, there is 
an equal amount of traffic generated destined for every other 
node and hence a limited amount of traffic destined for any 
particular destination. Under these assumptions we consider a 
particular message in the buffer of a particular node (call it the 
current node) and derive its expected waiting time for epidemic 
routing. Let p succ be the probability of successfully commu- 
nicating with the destination of the message given that the 
current node meets the destination (successful communication 
may not happen because of contention with other neighboring 
nodes). The number of meetings with the destination before 
the current node is successfully able to communicate with it 



is geometrically distributed with mean 



and the average 



time elapsed between meetings is Tj m . Hence, E[W] is given 
by 



E[W] = 



T 



(9) 



where p succ is given by 



where is the average number of communication sessions 
for which the current node remains in contact with the 
destination, and p succ is the probability that the current node 
successfully communicates with the destination in a particular 



communication session. p succ is given by 

oo 

Psucc = J2 P ( N »hb = iWnhb > 1) (PcomJ 2 

where p 3 comm is the probability that a particular node gets to 
communicate in a session given that it has j neighbors. p? comm 
depends upon the MAC protocol used to resolve contention. 
As described in Section IIII-A1 the MAC protocol effectively 
chooses all nodes in a cell if the number of nodes is atmost 
a, else it chooses a random a-sized subset of the nodes in the 
cell. Hence, 
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C. Contention Model 

We model contention in terms of three different factors as 
follows: (i) the probability of increase in number of copies of 
a particular message from i to i , namely p,-' ^ (it) the prob- 
ability of delivery to the destination of a particular message 
given that there are i copies of the message in the network, 
namely pd,i, and (Hi) the expected time lapse between the 
number of copies reaching i and either a rise in the number 
of copies or delivery to the destination, namely E[Di]. We 
assume that with a high probability, in any communication 
session, only at most one of the i nodes carrying a copy of 
the message attempts to transmit the message and hence pj i 
is zero for i > i + a. 

For i < i < i + a, since we assume that at most one node 
transmits this message in any communication session, 



a-l 

V P(N nhb = j) x min ( — r , 1 ] x 



j=i'—i 



3 + 1 

(fe«) i '- i x(l-p te (i)) i -( i '- i) 



where P(N n hb = j) has already been derived in Section IIV-AI 
and ptx{i) is the probability that a message is successfully 
transmitted to another node given they are in communication 
with each other and that there are i copies in the network. The 
probability ptx(i) is given by, 



Ptx(i) = 1 - 



i - 1 



/ b B w \ 



1J \E[B]J ■ 

To compute pd,i, we assume that at most one of the i nodes 
carrying a copy of the message meet the destination at a given 
point in time. Hence pd.i is the probability that any one of 
these nodes meets the destination and successfully transfers 
the message to the destination. And pd.i is given by, 



Pd,i 



3 = 1 



ntot - 1 



j + r 



a \ . / a 
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1) E[Di] Computation: To derive E[Di], we first estimate 
the average time elapsed before a particular node out of the 
i nodes carrying the message spreads a copy/copies of the 
message. We jointly model the process of arrival and departure 
of neighbors and the attempted transfer of the message as a 
stochastic process with states j = 0, 1, • ■ ■ n max representing 
the event that at the beginning of a communication session, 
the node has j neighbors and still has not spread copies of the 
message. The stochastic process also has a state C representing 
that the node has successfully spread one or more copies of 
the message and that therefore the number of message copies 
in the network has changed. This has been illustrated in Fig [2] 




n-system 

Figure 2: Digraph representing the stochastic process used to 
compute E[Di] 

We define the n-system as the subset of states n,n + 
1 • • • i nmax- Let E[t n ) be the expected dwell time in the state 
n, p s n n+ i(i) be the probability that the number of neighbors 
of the node under consideration increase from n to n + 1 
in a communication session with no message copies spread, 
given there are i copies of the message, n neighbors for that 
particular node, and that some change in the state of the 
stochastic process happens at the end of the communication 
session. Also, let p s s ^ n+1 ^ n (i) be the probability that a transi- 
tion happens from the (n+ l)-system to the state n given that 
a transition from the (n + l)-system occurs. The following 
lemma gives the expected dwell time E[t n ) of the process in 
the n-system. 



Lemma 4: 



E\t n ] = 



E[t n }+pi n+1 (i)xE[r 



X ^(n+l),„W 

Proof: We replace the dwell time distribution at each state 
n by its expected value E[t n ]. Then the expected time spent 
in the n-system is given by, 

=E[t n ]+ 

Pn,n+l(i) [Eitn+l] + P'l 

(n+l),n 

(i) x (E[t n ] + ■■■) 



E\f' 



which gives us 

E[t n ] = E[t n ] +K >n+1 (i) (E[t n+1 ] +p\ n+l)>n {i) x E[t n }) 

which can be rearranged to give us the result. ■ 
To terminate the recursion, we have 



To compute E[t n ] we first compute p s nn {i) which is the 
probability that the number of neighbors remains n at the end 
of a communication session with no copies of the message 
transferred, given that there are i copies of the message in the 
network. For n < a — 1 it is easy to derive 

Pl,„(i) = P{n\n) x pl(i) 

where P{n\n) is as derived in the aforementioned Sec- 
tion IIV-AI and is the probability that there is no suc- 
cessful communication of the message given that there are i 
copies of the message and n neighbors. For n < a — 1, p s n {i) 
is given by, 

pZ(i) = (i-p tx (i)) n . 

And for n > a — 1, p s n {i) is given by, 

Pi(i) = (i--+-x(i- Ptx (i)r). 

V n n / 

The number of communication sessions for which the system 
stays in state n is then a geometric random variable with 
probability of success 1 — p s n n (i). Hence, the expected dwell 
time in state n is given by, 



E[tn] = - 



Also, it must be noted that since the session time is small, 
the expected dwell times E[t n ] can be approximated by 
continuous exponential distributions. Hence, if we considered 
i nodes in state n instead of one and computed the expected 
dwell time (the dwell time in this case would be the minimum 
of the dwell times of the i nodes), the expected dwell time 
would approximately be 



E\t„ 



The transition probability p s n n+1 (i) is given by, 



P 



n,n+ 



p* n (i)x P(n + l\n) 
1 - pf^i) x P(n|n) ' 



To compute the other transition probability p*/ n+1 N n (i), we 
assume that the dwell time within the (n + 1) -system is 
sufficient for the stationary distribution within the (n + 1)- 
system to be reached. Hence we have, 



P S s (n+l),nW 



v n +i x P(n\n + 1) 



x P(n\n + 1) + Ek=n +1 v k x(l- p°(i)) ■ 



The expected time elapsed before a particular node carrying 
the message spreads copies of the message is finally given by 
E[t°}. However, E[Di] is given by the expected time elapsed 
before any one of the i nodes carrying the message spreads 
copies of the message. Since each E[t„] is approximately 
inversely proportional to the number of nodes we consider 
for the stochastic process, and E[t°] is a linear combination 
of E[t n ], we can approximately say 



E\t n " 



E[t, 



E[Di 



E[t°] 



V. Routing Delay 

In this section we derive the routing delay for epidemic 
routing in a DTN which uses M2M communication. We model 
the evolution of the network with respect to a specific message, 
from the generation of the message to the spreading of copies 
of the message to finally the delivery of the message, as a 
stochastic process, represented as a digraph in Fig [3] This 
stochastic process has a set of states i = 1, • • • , n to t — 1 where 
state i corresponds to the presence of i copies of the message 
in the network, and ntot refers to the total number of nodes 
in the network. The stochastic process also has an absorbing 
state D which represents successful delivery of the message 
to the destination. The network spends an average time of 
E[Di] in state i before moving on to one of the next possible 
states. Since we assume no message drops (i.e., buffers are 
adequately sized), the probability of transition p,i i from state 
i to i is zero for i < i. When the network is in state i, we 




Figure 3: Digraph representing the stochastic process used to 
compute the routing delay 



assume that the probability that more than one of the nodes 
carrying a particular message try to transmit the message at 
the same time is negligible. Hence, when the network is in 
state i, at most one of the i nodes will transmit the message. 
This can be transmitted to at most a — 1 other nodes (since 
at most a nodes can participate in an M2M communication 
session). Hence, the network can evolve from state i only to 
states i + 1, • • ■ , i + a — 1 . We call the set of states reachable in 
a single transition from state i as Ch(i). This corresponds to 
the children of node i in the digraph representing the stochastic 
process (assuming edges corresponding to zero-probability 
transitions are deleted). Similarly, a state i can be entered only 
from states i — 1, • • • , i — a + 1. We call the set of states from 
which state i is reachable in a single transition as Pa(i). This 
corresponds to the parents of node i in the digraph. 

We then derive the expected delay E[D] of epidemic routing 
by averaging the delay over all the possible paths $ on the 
digraph starting from the starting state 1 and ending in state 
D. It must be noted such a path corresponds to a sequence of 
states si, s 2 , • • • , s m where Si = 1 and s m = D. 

E[D] = ]TP ($)£($) 

= Yl (P(^\l)---P(D\s m ^))x 
(E[D 1 ] + ---E[D m _ 1 \). 



Table I: Simulation settings 



Parameters 


Values 


Terrain size (in mxm) 


100 x 100 


Grid size (in m) 


10 x 10 


M2M communication session time (in seconds) 


0.1 


M2M a-value 


4 


Message size (in Kb) 


1 


Message generation rate 


1 

5 


Message TTL 


CO 


Buffer size 


CO 


Scheduling 


FCFS 


Simulation time (in seconds) 


10000 



This expression can be re-written as follows: 

^] = E E E 

s pePa(s) ceC'h(s)U{D} 

P(Path(l,p)). Ps , p .p c , s .P(Path{c,D)).E[D s ] 

where P(Path(i, j)) represents the total probability of all 
paths starting from state i and ending at state j, given that 
the system is initially in state i. P(Path(l, i)) is given by the 
following recursion 

P(Path(l,i)) = P(Path(l,p)). Php . 

pEPa(i) 

And this can be calculated for all states i efficiently using 
dynamic programming. Similarly, P(Path(i, D)) is given by 
the recursion 

P(Path(i,D)) = J2 {pc,i-P(Path{c,D)))+p Dti . 

cech(i) 

And this can also be calculated efficiently for all i using 
dynamic programming. 

VI. Performance Evaluation 

The simulations were done using the ONE simulator |[T3l . 
which we custom modified to simulate one-to-one CSMA 
and M2M communication schemes. The simulation settings 
are shown in Table Q] Each experiment was run for 10000 
simulation seconds. The system was verified to have reached 
the stability region within this time. As mentioned earlier, 
traffic generation is symmetric, i.e., every message generated 
is given a randomly chosen destination node. 

Figure [4] compares the routing delay values obtained from 
analysis and simulations for 10 and 20 nodes, and for two 
different link bandwidths, namely 15 kbps and 50 kbps. For 
each configuration, the comparison is done over average node 
speeds (zT) ranging from 5 m/sec to 8 m/sec. These delay values 
were observed in the stability region of the experiments and 
the delivery ratios (i.e., the fraction of messages delivered 
successfully) were observed to be close to 0.9. While we have 
made a number of approximations in developing our analytical 
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Figure 5: Routing performance of CSMA and M2M communication schemes for 100 nodes with two different average node 
speeds (v), namely 7 and 8 m/sec 
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Figure 4: Delivery delay performance over two different 
number of nodes (namely, 10 and 20) and across different 
average speeds (v), for two different link bandwidths 



model, we find that the delay values from analysis and 
simulations match, thereby validating our model and showing 
that the approximations were reasonable ones to make. As is 
expected, it is observed that with an increase in node speeds 
the routing delay decreases, as the mixing of nodes enable 
faster delivery of messages. Moreover, an approximately three- 
fold increase in the link bandwidth proportionally reduces the 
routing delay, which is intuitive as messages are pumped three 
times faster for the same contact time scenarios. 

Figure [5] compares the routing performance of one-to-one 
CSMA and M2M schemes. These experiments were run with 
a lower link bandwidth of 10 kbps for M2M and 40 kbps for 
one-to-one CSMA. CSMA is given a higher bandwidth since it 
uses the entire spectrum available while M2M communication 
divides the spectrum among the a (= 4) nodes participating 
in a communication session using FDMA. The M2M scheme 
clearly outperforms the one-to-one CSMA scheme. We find 
that when the M2M scheme reaches its stability region, the 
CSMA scheme has not yet attained stability, and the differ- 
ences between the two schemes are increasing with time. The 
fact that M2M has reached its stability region while CSMA has 
not, is supported by the high delivery ratio values (i.e., above 
0.9) for M2M and lower delivery ratio values (i.e., around 
0.8) for CSMA. Also our M2M delay values stabilize near 
the corresponding analytical delay values. 

VII. Conclusion 

In this paper, we have shown for the first time how M2M 
communication can be used to combat contention, and thereby 
reduce the routing delay in a DTN. A theoretical model has 
been developed for the delivery delay of a DTN employing 
epidemic routing and M2M communication and has been 



validated against simulations. Using simulations, we have also 
shown that M2M communication significantly outperforms 
communication that uses one-to-one CSMA in terms of deliv- 
ery delay. In the future, we plan to extend the model developed 
to work with more sophisticated mobility models. 
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